How well can (renormalized) perturbation theory predict dark matter clustering 

properties? 
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There has been some recent activity in trying to understand the dark matter clustering properties 
in the quasilinear regime, through resummation of perturbative terms, otherwise known as the 
renormalized perturbation theory or the renormalization group method 2J. While it is not 
always clear why such methods should work so well, there is no reason for them to capture non- 
perturbative events such as shell-crossing. In order to estimate the magnitude of non-perturbative 
effects, we introduce a (hypothetical) model of sticky dark matter, which only differs from collisionless 
dark matter in the shell-crossing regime. This enables us to show that the level of non-perturbative 
effects in the dark matter power spectrum at k ~ 0.1 Mpc -1 , which is relevant for baryonic acoustic 
oscillations, is about a percent, but rises to order unity at k ~ 1 Mpc . 



In the era of precision cosmology, an accurate under- 
standing of dark matter clustering properties is an es- 
sential ingredient of relating observations of clustering in 
galaxy or weak lensing maps to the fundamental prop- 
erties of our Universe, such as its geometry or linear 
growth history. This need is underlined by the recent 
discovery of baryonic acoustic oscillations in the power 
spectrum of galaxies [3, 0] , and the prospects of its ap- 
plication as a cosmic standard ruler (e.g. ||). While 
semi-analytic methods, most famously within the con- 
text of the halo model [f| 0, Q , have been successfully 
used to model both galaxy and dark matter clustering 
properties, they suffer from their phenomenological na- 
ture, which inevitably leads to a plethora of parameters 
that need to be calibrated using numerical simulations or 
observations. Moreover, the number of these parameters 
is expected to increase through more in-depth studies, 
which will certainly lead to the surfacing of more sub- 
tle effects (e.g. g QjJ). Thus, the ultimate dilemma 
in the era of precision cosmology will become if we are 
probing the fundamental properties of our universe, or 
rather further constraining the intricate phenomenology 
of non-linear gravitational gastrophysics. 

Numerical simulations go a long way in elucidating the 
complicated process of gravitational collapse of dark mat- 
ter and baryonic gas. However, they still suffer from our 
poor understanding of star and galaxy formation, and 
its feedback on the surrounding intergalactic medium. 
Moreover, they may also suffer from finite resolution and 
box size effects, as well as potential ill-understood nu- 
merical artifacts that can limit the accuracy of numerical 
studies. 

The latter has motivated analytic studies, which are 
most systematically done in the context of perturbation 
theory (see [ll| for a review). However, consecutive terms 
in a perturbative series become comparable in the quasi- 
linear regime, when overdensities become of order unity, 
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signaling the breakdown of the perturbative expansion. 
Inspired by the use of renormalization methods in quan- 
tum field theory, some studies [l|, 0] have introduced 
renormalization or resummation technics that, in prin- 
ciple, capture the dominant terms in the perturbative 
series, and can yield an accurate result, even when the 
perturbation theory breaks down. 

While the theoretical validity of the technics advocated 
in these studies is not transparent to this author, it is 
possible to put a limit on the degree of accuracy that any 
method, based on perturbation theory, can achieve. This 
limit comes from the non-perturbative effects involved in 
the gravitational collapse process, and is the subject of 
this note. 

The cosmological perturbation theory for cold dark 
matter [ll[ is based on the pressureless Euler equation, 
where the fluid starts from a nearly homogenous expand- 
ing initial conditions, and then evolves under its own 
gravity. The equations break down at shell crossing, 
which is when perturbations grow to the limit that mul- 
tiple streams cross each other, and thus the assumption 
of zero pressure fails. 

As long as different streams of dark matter do not cross 
each other, the cross-section for their collision does not 
enter the equations. Therefore the perturbative frame- 
work, which is the basis for the recently suggested renor- 
malization methods, is insensitive to dark matter colli- 
sional properties. Even though in the minimal dark mat- 
ter model, the particles have a negligible self-interaction 
cross-section, one may also envisage an opposite regime, 
where the collisions have a huge cross-section, and are 
maximally inelastic. We dub this particular model for 
dark matter as sticky dark matter. 

Of course, sticky dark matter is a terrible candidate 
for the cosmological dark matter. This is because, as 
a result of frequent and inelastic collisions, all the dark 
matter particles would sink into the center of a halo, 
probably forming a black hole, or a small compact disk (if 
prevented by the initial angular momentum of the halo) . 
However, it provides an interesting theoretical test study 
to examine the non-perturbative effect of shell-crossing, 
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FIG. 1: Left: Halo model dark matter power spectrum for collisionless (solid) and sticky (dotted) dark matter. Right: Relative 
difference between power spectra of collisionless and sticky dark matter models. This shows the level of non-perturbative effects 
in the dark matter power spectrum. 



as sticky and collisionless dark matter obey the exact 
same equations within the (single-stream) perturbation 
theory. 

In order to make this comparison, we resort to the 
simplest version of the halo model , where dark mat- 
ter haloes cluster according to the linear matter power 
spectrum with a constant bias, which only depends on 
their mass. The structure within a halo of collisionless 
dark matter is well approximated by an NFW form [12j: 

P{r) = (r/r 8 )(l + r/r s )*' W 

which extends out to r„j r = cr s , where the concentration 
parameter, c, as a function of the total enclosed mass is 
also constrained from simulations. 

It is known that a combination of halo-halo corre- 
lations, plus single halo auto-correlation terms gives a 
reasonable approximation to the dark matter non-linear 
power spectrum in numerical simulations (e.g. [lH). 
With this qualitative picture at hand, we can thus ex- 
amine the case of sticky dark matter, by collapsing each 
NFW halo into a point. 

Fig. ([T]) compares the dimenionless power spectra, 

A 2 (fc) = k ^Jii^ i f° r collisionless and sticky dark mat- 
ter models, in the context of the halo model. Here, we 
assume the WMAP 1 st y ear best fit concordance cos- 
mological parameters [14]]. We see that the difference 
between the two models, which characterizes the level of 
non-perturbative effects in the dark matter power spec- 
trum is ~ 1% at k ~ 0.1 Mpc" 1 , which is the region 
relevant for probing baryonic acoustic oscillations. This 
level rises to order unity at k ~ 1 Mpc -1 , indicating the 



complete breakdown of (renormalized) perturbation the- 
ory, or any single-stream renormalization group method. 

Fig. ([2]) shows the same difference between the sticky 
and collisionless power spectra, as a function of the 
non-linear power. It is interesting to notice that non- 
perturbative effects do not dominate the power until 
A 2 (fc) ~ 100, which is far into the non-linear regime. 
This may explain why renormalization methods based 
on the perturbation theory, such as [l|, [2J , can do so well 
in the quasilinear regime. 

A question that may arise is if we are stretching the 
halo model beyond its level of applicability. It is true that 
the halo model has had significant phenomenological suc- 
cess in fitting simulations and observations (e.g. [l3j and 
[To]]). However, it suffers from unphysical features, such 
as spurious power on very large scales, which is caused 
by the presence of the 1-halo term @, [3]- However, the 
difference shown in Figs. ([I]) and ((2|) is dominated by the 
2-halo term for k < 0.3 Mpc -1 (at least down to near- 
horizon scales). Even without any reference to the halo 
model, it is clear that the large scale correlation function 
is, to some extent, affected by the change in the profile 
of the virialized regions, say from an NFW (for collision- 
less dark matter) to a point-like profile (for sticky dark 
matter). The difference in the 2-halo term, in the con- 
text of the halo model, gives an approximation for the 
magnitude of this effect. 

As vorticity cannot be locally generated in the single 
stream regime of dark matter collapse, a different test 
for the breakdown of the single-stream perturbation the- 
ory is when the vorticity and the divergence of the ve- 
locity field become comparable [l?) ■ This happens at 



3 




1Q -5 I I 

0.001 0.01 0.1 1 10 100 

A 2 (k) 

FIG. 2: The level of non-perturbative effects in the non-linear 
power spectrum (see FigfT} , as a function of non-linear power. 

k ~ 2 Mpc -1 , which is consistent with our results. How- 
ever, in contrast to our approach, it is not clear how 
the vorticity to divergence ratio is related to the level 
of non-perturbative effects in the power spectrum within 
the quasilinear regime. 

What if a single-stream perturbative method agrees 
with numerical N-body simulations at a better level than 
predicted by our simple exercise? Then I would argue 
that this level of accuracy is not warranted for any single- 



stream calculation, and thus any such agreement is acci- 
dental, unless it is justified based on a non-perturbative 
method. 

Finally, we should point out that sticky dark matter is 
very similar to the so-called adhesion model [18|, [l9( , that 
has been introduced in the context of Zel'dovich approx- 
imation [2(3]. Adhesion model, which is realized as the 
zero/small viscosity limit of a (single stream) fluid, has 
been introduced to enable the use of analytic Zel'dovich 
approximation beyond shell-crossing (or pancake forma- 
tion). However, direct comparison of the adhesion model 
to the sticky dark matter might be misleading, as it mixes 
the perturbative terms, missing in the Zel'dovich approx- 
imation, with the non-perturbative adhesive feature. 

In conclusion, in this note, we have pointed out that 
any renormalized perturbation theory method is ulti- 
mately limited to the context of single stream fluid equa- 
tions, and thus does not capture the non-perturbative 
shell-crossing events. To estimate the level of non- 
perturbative effects, we introduced a toy model of sticky 
dark matter, which only differs from collisionlcss dark 
matter in the shell-crossing regime. This provides us 
with a framework to estimate the level of accuracy ex- 
pected from renormalized perturbation theory, showing 
that non-perturbative effects could be ~ 1% at scale of 
baryonic acoustic oscillations, but do not dominate until 
deep in the non-linear regime. 
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